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Abstract. - 

We present generalized Green-Kubo expressions for thermal transport coefficients /i in non- 
conservative fluid-type systems, of the generic form, fi — /ioc + J °° dtV- 1 (I e exp(t£)/) where 
exp(i£) is a pseudo-streaming operator. It consists of a sum of an instantaneous transport 
coefficient fioo, and a time integral over a time correlation function in a state of thermal 
equilibrium between a current / and its conjugate current 7 e . This formula with /j.^ ^ 
and covers vastly different systems, such as strongly repulsive elastic interactions in hard 

sphere fluids, weakly interacting Langevin fluids with dissipative and stochastic interactions 
satisfying detailed balance conditions, and "the likes", defined in the text. For conservative 
systems the results reduce to the standard formulas. 



The Green-Kubo formulas for thermal transport coefficients in simple classical fluids with 
smooth conservative interactions are widely used, and generally accepted [1, 2, 3] as exact 
expressions for general densities, as long as the deviations from equilibrium and the gradients 
are small, and the transport coefficients exist. This standard expression is given in terms of 
an equilibrium time correlation function between iV-particle currents I and l t , i.e. 



where fi denotes a typical transport coefficient, (. . .)q is an average over a thermal equilibrium 
ensemble at temperature T = l/fee/3, with a distribution function po ~ exp[— /3H], and 
limv^oo denotes the thermodynamic limit. In the sequel this limit is understood, but not 
explicitly written. For convenience, time is taken to be continuous. The time evolution 
of a dynamical variable, I(t) = e tc I(0), can be described by a streaming operator e tc , 
which generates a time evolution that is invariant under the time reversal transformation. 
Consequently, the Liouville operator satisfies, £ c = where C e denotes the time reversal 
transform of C. As C and the currents / contain in general interparticle forces, the standard 
Green-Kubo formula (1) is ill-defined if the interaction potentials are not sufficiently smooth, 
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and it does neither apply to hard sphere or hard core interaction potentials (non-conservative 
impulsive forces), nor to dissipative and stochastic forces, which are not derivable from an 
interaction potential. 

Problems addressed: The goal of this letter is to present generalized Green-Kubo expressions, 
without any restrictions to dilute systems, for a large class of complex (non-conservative) fluids, 
i.e. -/V-particle systems, that do approach a state of thermal equilibrium, but whose equations 
of motion have ill-defined forces, or/and lack time reversal invariance. The forces may be 
impulsive, dissipative, stochastic, and possibly include conservative forces as well. The generic 
form of this generalized Green-Kubo formula is, 



where it must be possible to represent the time evolution through a pseudo-streaming operator 
exp(tC), which is well-defined for all points in the relevant phase space, and generates the 
proper trajectories of the dynamical variables when used inside statistical averages (• • -) . The 
prefix pseudo refers to the fact that the representation exp[i£] is only valid inside averages. 
An assessment of the type of systems to which this Green-Kubo formula is also applicable, 
called " the likes" , and an outline of the derivation and implications are the main goals of the 
present article. The discussion covers such vastly different systems as microscopic fluid models 
with strongly repulsive elastic interactions (hard spheres, square well potentials) for which the 
interparticle pseudo- forces and the pseudo-Liouville operators in Eq.(2) are well-defined, as 
well as mesoscopic weakly interacting Langevin fluids with dissipative and stochastic forces. 

There seem to exist in the literature only three cases where explicit Green-Kubo formulas 
of the generalized form (2) have been derived for transport coefficients in complex fluids, each 
referring to a very different system and to a very different transport coefficient, obtained by 
very different approaches. The first case in Ref.[4] refers to Brownian dynamics. Here the 
friction coefficient of a Brownian particle in a fluid, both modeled by hard spheres, is derived 
with the help of a delicate analysis, using an asymptotic expansion in the small ratio m/M 
of the hard sphere mass over the mass of the Brownian particle. In this limit the interactions 
between the heavy particle and the liquid particles becomes weak, and the kinetic equation 
for the distribution function of the heavy particle reduces to a Fokker Planck equation. Here 
the friction coefficient is expressed as a generalized Green-Kubo formula (2), and contains a 
time correlation function between a pseudo-interparticle force and its conjugate. A possible 
extension of this method to complex fluids is not obvious. 

A second case, also concerning hard spheres, can be found in Rcf. [5] for the shear viscosity 
of an elastic hard sphere fluid, derived by a transformation of the corresponding Einstein- 
Helfand formula, or by a limiting procedure applied to the standard Green-Kubo formulas for 
soft spheres with a repulsive interaction potential V(r) ~ 1/r" with n — > oo. The limiting 
procedure cannot be extended to complex fluids, but the first route via an Einstein-Helfand 
formula is possible, following Ref. [5]. 

A third case concerns a Langevin fluid [6] , where standard linear response theory [1] is used 
to obtain the long time diffusion coefficient of a colloidal suspension in the generic form (2). 
This is done by modeling the short time diffusion coefficient by a phenomenological relaxation 
constant, used as input in the Langevin equation. The extension of that method for a weakly 
interacting Langevin fluid to strongly interacting hard sphere systems does not seem feasible 
either. 

The question of interest in this article is then: do these scattered results have something 
in common that makes the structure (2) generic for a larger class of physical systems? The 
answer to this question will be given below. 
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Properties of £: Suppose that £ has been constructed. As we are studying equilibrium time 
correlation functions, we have to impose the condition of stationarity, 

(A(0)B(t)) = (A(-t)B(Q)) . (3) 

Here the time evolution of A(— t) = cxp(t£ e )A(0) (where t > 0) is generated by the backward or 
time reversed streaming operator. Defining the transpose £ through the relation, J dTAe tc B = 

J dTBe tc A, we find by comparing integrands in Eq. (3) the necessary condition on £ in the 
form of a commutation relation, i.e. 

p Q £ = £ e p or Z = Z\ (4) 

where £ is defined through po£ = £po- The above ingredients arc sufficient to derive 
Green-Kubo type formulas for the Navier-Stokes transport coefficients using linear response 
theory. The class of "the likes" consists of those microscopic or mesoscopic models that have: 
(i) a pseudo-streaming operator of the form exp[i£], (ii) a stationary distribution given by 
the thermal po, where (iii) pa and £ satisfy the commutation relation (4). Note that for 
conservative interactions, where £ = —£ = £ e 7 the condition (4) is trivially satisfied as 
£ = £ £ — £, and £ and po commute. 

There exist two important simple realizations of the previous scenario. The first and most 
important class are the fully microscopic models of hard sphere fluids, which are prototypical 
for the strongly repulsive hard core interactions of classical fluids. The other class are meso- 
scopic Langevin fluids with weak, dissipative and stochastic forces, where the interactions are 
modeled by phcnomcnological input parameters. 

For hard spheres the standard Liouville operator for conservative systems is ill-defined 
because of the impulsive interaction forces. For this case the construction of pseudo-Liouville 
operators £ in terms of binary collision operators is far from trivial [7]. The resulting 
expressions for £ in terms of binary collision operators satisfy the commutation relation (4), 
as shown in [7]. This pseudo-Liouville operator is not odd in the velocities. So £ e ^ —£, and 
consequently not all trajectories generated by exp(i£) are time reversal invariant. 

For Langevin fluids the construction of the pseudo-Liouville operator £ is rather trivial, 
as this ^operator can be simply obtained from the corresponding Fokker-Planck equation, 
dtp = £p, for the TV-particle probability distribution p(T,t) [8]. As this equation satisfies 
the detailed balance condition, the stationary solution of dtpo = £po = is the thermal 
distribution po. Also here it is sufficient to define pseudo-streaming operators e tc for the time 
evolution of dynamical variables only inside averages. This may be done through the relation, 

(A(t)) = J dTAe tC p{Y, 0) = J dT P (T,0)e tc A=(e t£ A). (5) 

Here the pseudo-Liouville operator £ is the transpose of the Fokker-Planck operator. In this 
formulation £ is not the infinitesimal generator of the Langevin equations, but an effective 
operator, that acts inside averages, in which the rapid fluctuations of the random forces have 
been averaged out. Because of the presence of dissipative forces, the Fokker-Planck equation 
is not invariant under the transformation of time reversal. Consequently, £ e ^ —C. For 
Fokker-Planck equations, satisfying the detailed balance criteria, the operator £ satisfies the 
commutation relation (4) (sec [8]). 

Linear response theory: Having discussed the properties of the pseudo-Liouville operator, 
one can apply linear response theory for calculating thermal transport coefficients in fluids 
[1, 3]. Consider, for simplicity, a system with a single conserved density e(r). Extensions 
to more conservation laws are straightforward. Let ek be the Fourier mode of the conserved 
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density, ek = i5ejexp[— ik • r*] with fluctuation <5e, = e — (e} . The microscopic equation of 
motion for the Laplace transform ekz of ek(i) is (z — £)ek2 = ek- In the long wave length limit 
(k — > 0) the hydrodynamic propagator satisfies the relation, 

G(k, z) = (e k |e kz ) ~ (eo|e ) /[z + k 2 D], (6) 

where D = /j,/ (e |e ) is the diffusivity, and \i the kinetic coefficient. The standard projection 
operator method [1] enables us to derive in a few lines an exact expression for D or \i. To do 
so we introduce the projection operator V = 1 — V±_ acting on a dynamic variable 6k as Vb^ = 
ek(ek|frk)/(ek|ek)- The inner product is defined as a thermal average, (ak|& k ) = V^ 1 (a£&k) - 
Application of this method yields, 

[z - VLV - VCVAz - V ± CV ± )- 1 V ± CV]Ve kz = e k . (7) 

Writing this in component form, and comparing the result with (6) yields for D, 

D = -Cim ft [(e k |£e k ) + (£ £ e k |7>j_&7>]£e k )] /[k 2 (e k |e k >], (8) 

where Cim represents the double limit, lim z ^ lim^Oj denotes the real part, and Q z reduces 
in the small-Zc limit to the projected resolvent, V±(z — C)~ l T'_\_. To obtain the second term we 
have used the relation (e^\CV± ■ ■ •) = (C t ek\P± • • •}, based on the commutation relation (4). 
This term can be expressed as a current-current correlation function. 

Consider the so called Euler matrix (e k |£ek) = — ik (ek|i k ). Its imaginary part would yield 
the Euler equations for a fluid with the standard conservation laws. As conserved densities 
have in general a definite parity in the velocities (mass, momentum, energy), the Euler matrix 
element would vanish for conservative systems, where C e = —£. However, for systems whose 
pseudo-streaming operator lacks time reversal invariance, C e ^ —C, and the matrix element 
has a non- vanishing real part for small k of 0(k 2 ), i.e. 

Moo = (e |e ) Doo = - lim fc~ 2 5ft (e k |£e k ) . (9) 

k — >oo 

Next consider the second term in Eq. (8). Here the total microscopic flux, J = lim^o j k , is 
related to the Fourier mode ek of the conserved density through the local conservation law, 
c*te k = £e k = —ikjk. Similarly, we define the conjugate current, J e = lim.k-yo j*u through 
£ e a£ = —ikj* k . Inserting these definitions into Eq. (8) yields finally for the kinetic coefficient 
\x the generalized Green-Kubo formula (2), where / and I e are subtracted currents, defined as 
I = V±J. This completes the formal derivation of Eq. (2), and it remains to work out the 
explicit expressions for /, I e and ^oo. In the remaining part of this letter we present some 
applications to a heat conducting random solid, an isothermal Langevin fluid, and a hard 
sphere fluid. 

Heat conducting random solid: A realization of a system in the previous section is a quenched 
Langevin fluid with heat conduction. It consists of N point particles, quenched in a random 
configuration X = {r^ji = 1, • • • , N}. Each particle is characterized by its internal energy ej, 
collectively denoted by e = {ei{t)\i — 1, • • • , N}, which are the only dynamical variables in the 
model. The internal energy corresponds to many internal degrees of freedom, described by a 
density of states ~ e a , where a is proportional to the number of internal degrees of freedom 
(a 3> 1). Energy is exchanged between the particles through dissipative and stochastic forces of 
finite range, which obey detailed balance criteria. The total energy, E = ^ e,(t) is conserved, 
and the A^-particle system approaches a state of thermal equilibrium. The system supports a 
local heat current, driven by gradients in the local energy density or temperature field. So, at 
the macroscopic level Fourier's law of heat conduction applies, at least if the density of particles 
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is above a critical percolation threshold. Our goal is to derive a Green-Kubo expression for its 
heat conductivity. 

The Langevin equations for the time evolution of the dynamical variables ej(t) and the 
corresponding Fokkcr-Planck equation for the TV-particle probability density p(e, t\X) in a 
fixed configuration X of N point particles have been derived to dominant order 0(1/ a) in 
Ref. [9, 10]. We only quote the adjoint Fokker-Planck operator £ which reads, 

c = T,i<j x ij [i e 3 - e i) + a^eiejdij] dij = C e . (10) 

Here dij — djdti — d/dej 7 and \j = X w(rij) with a positive range function w(rij) vanishing 
for rij > r c , and 70 is a phenomenological model parameter. The equation above shows that 
C e — C because €i does not change sign under a time reversal transformation. The microscopic 
energy fluxes for small k are obtained as [10], 

£e k ~ ik ■ J2i<j^ij( e i ~ e i) r i] = ' Q 

C e e* k ~ ik ■ Q £ = Cel = +ik ■ Q, (11) 

and there are no subtracted parts. Note the relation Q £ = — Q. The Green-Kubo formula for 
heat conductivity in the random solid is then A = Aoo + X dd , with 

/>oo /*oo 

\ dd = [fil dTV) I dt(Q,-e tc Q) = -(P/dTV) dt(Q-e tC Q) Q . (12) 
Jo Jo 

The Green-Kubo expression (12) has the generic form for systems with dynamics lacking time 
reversal invariance. As this formula can also be derived from an Einstcin-Hclfand formula, 
A is necessarily positive. Because £ £ ^ —£, necessarily 7^ 0. The instantaneous heat 
conductivity is equal to the mean field approximation to the heat conductivity, calculated in 
Refs.[9, 10], i.e. A m = (\on 2 C v /2d)[wR 2 ] where [wR 2 ] = j (iRw(R)i? 2 and C v is the specific 
heat per particle. In fact, the remaining time integral in (12) gives at large densities only a 
negligible contribution, but at the percolation threshold both terms cancel [10]. If the standard 
Green-Kubo formula (1) would have applied, then A M = 0, and the sign in front of the second 
time integral would have been a plus sign. So, the resulting expression would have been quite 
different from Eq. (12). 

Isothermal DPD fluid: Next we will consider the so called isothermal DPD (dissipative 
particle dynamics) fluid, introduced in Rcf.[ll], and studied analytically [12, 13] and by 
means of computer simulations [14]. It is a mesoscopic version of a classical fluid where 
the fast microscopic length and time scales are averaged out, and DPD particles, described 
by their position rj and velocity , can be thought of as soft lumps of fluid. Their equations 
of motion, dvi/dt = Vj and mdwi/dt — ^^(F^ + + F^), contain conservative forces, 
F?j = —dV(rij)/drij, as well as dissipative (d) and random (r) forces, where F d is proportional 
to a phenomenological relaxation parameter 70. Explicit expressions for F d and F r can be 
found in Refs. [12, 13]. All forces have a finite mesoscopic range r c . The corresponding 
Fokker-Planck equation satisfies the detailed balance condition. The pseudo-Liouville operator, 
C = C c + d, is obtained as the transpose of the Fokker-Planck operator. Its conservative part 
C c includes the inertial term and the conservative forces (where C e c = —C c ). Its dissipative 
part Cd contains the combined action of the dissipative and stochastic forces (where C d = Cd)- 
Consequently the instantaneous transport coefficient depends solely on the dissipative force. 
In the DPD fluid total mass and momentum are conserved, but not the total energy. So 
mass density nk and momentum density gk are the only slow Fourier modes. Therefore, the 
fluid exhibits shear and bulk viscosity. The transport coefficients are given by the matrix 
generalization of Eq.(8) (see [3]). Here we restrict ourselves to the shear viscosity 77. It is 
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given by an expression similar to (8) with ek replaced by the transverse momentum density 
<7k_i_ = J2i mv iy ex P[ — ikr ix \ where k is taken parallel to the x-axis for convenience, and it is 
straight forward to obtain, 

/>oo 

v = voo + ip/v) / dt((r xy J d xy )e tc {.r xy + j d y )) o 

Jo 

= Voo + V cc + V cd + V dc + V M (13) 

For small k the relation £<7k_L — —ikJ xy has been used, where the total momentum cur- 
rent J xy = J xy + J xy , and J exy — J xy — J xy . Here J xy is the momentum current in the 
standard Green-Kubo formula for conservative forces, given here by t] cc , and J xy = — J d xy = 
Si<j r ij.xF d j y is the dissipative current. The instantaneous viscosity ij^, which is proportional 
to fc~ 2 (gk±\Cgu±), is equal to the quantity r\D = mn 2 -fo[gowR 2 ]/2d(d + 2), where go{r) is the 
radial distribution function in thermal equilibrium, calculated in Ref.[ll, 13] as the mean field 
approximation to the viscosity. Also note that the commutation relation (4) implies that 
rf d — i] dc 7 which contributions are in general non- vanishing. If the standard Green-Kubo 
formula would apply, then rj^ = 0, and the minus sign in the integrand would be a plus sign. 

Hard sphere fluid: We start from the expressions for C = L + and C e = —L- in the notation 
of hard sphere kinetic theory [7, 15, 16, 5], i.e. 

L± = Lo± E^.TMij); L ± - L ± E»<^±fe')- ( 14 ) 

Here L = Ej v i - d/dri is the inertial part. The binary collision operators satisfy the relations 
T± = T± and T± — T T , and their explicit forms can be found in the literature [7]. Then the 
linear response theory leads to the generalized Green-Kubo formulas (2) for the Navier-Stokes 
transport coefficients in hard sphere fluids, which reads e.g. for the shear viscosity, 

V = Voo + (0/V) f dt (j_ xy e tL +J +xy ) Q . (15) 
Jo 

Here J xy is the total momentum current, obtained for small k from the relation L±<?k_L — 
—ikJ± xy = —ik(J xy + J± xy ), where (fc) refers to the kinetic flux, J xy = ^ Ji mVi X Vi y , and (v) 
to the collisional transfer flux, 

J ±xy = ± Y.i<j T ±( i j) m { v iyn x +v 3V r JX ) 

= Ej<j mcr<i / (=F) d ^( v ij ■ o fd x o y 5{v lj - aa). (16) 

The constraint (=p) on the a-integration restricts that integral to the pre/post-collision hemi- 
sphere where TVy ■ a > 0. The explicit expressions for the currents, which are identical to 
those in Ref.[5], are very different from the ones appearing in the standard form (1), and 
also very different from those in those in the Langevin fluids. The instantaneous viscosity 
is proportional to k~ 2 (<jrk_i_|£s , k_i_) , and is calculated in Rcf. [5] as = pcr 2 /[d(d + 2)t^], 
where a is the diameter of the hard spheres, and the Enskog mean free time. The same 
instantaneous contributions can also be identified [15, 17] in the Enskog theory [18] for the 
hard sphere fluid. Details are presented in [17]. 

We conclude with a number of comments: 
(i) We have derived generalized Green-Kubo formulas, Eq. (2), for non-conservative systems, 
ranging from microscopic hard sphere fluids with elastic impulsive interactions, to mesoscopic 
Langevin fluids with weak dissipative and stochastic interactions. Moreover a fundamental 
commutation relation (4) has been formulated which guarantees that the time correlation 
functions are stationary, and characterizes those systems for which the transport coefficients 
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are given by the generalized form (2). For conservative systems the generalized formulas reduce 
to the standard Green-Kubo expressions (1). 

(ii) The generalized formula (2) does not apply to fluids of inelastic hard spheres, which 
are models for granular fluids, although infinitesimal generators of the form (14) have been 
constructed. The reason is that such dissipative systems, cither freely cooling or driven by an 
energy source, do not reach a stationary state described by the thermal distribution p . 

(iii) From the existence of pseudo-streaming operators for hard sphere fluids an interesting 
side result can be deduced. It allows us to introduce pseudo-interparticle forces [4, 5], F^j 5 
and conjugate forces through the relation CmVi = J2j(^i) ^f/i an d £ e mvi = J2j(^i) ^e,ij 
respectively. An identity, derived in [5], states that (F^- x (0)F^j y (t)) = CS xy S(t), where 
C is given explicitly. Comparison with the time correlation function for the Langevin forces 
suggests that the pseudo-interparticle forces in hard sphere systems can be considered as white 
noise. 

(iv) Another side result of interest, following from (14), is the observation that the frequently 
cited Grccn-Kubo formulas for the DPD fluid [19], regarding for the viscosity for instance, 
r\ = i] cc — rf d , arc not correct (compare with (13)). Equally incorrect is the standard 
Green-Kubo formula (1) for the heat conductivity of the random solid, quoted in Ref.[10] 
as A = — \ dd (compare with (12)). It was mentioned there as a possible alternative for the 
kinetic theory analysis, presented in that paper. Regarding the main interest of this letter, 
i.e. the generic structure of the Green-Kubo formulas for non-conservative systems, the last 
comment is only of minor importance. 
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